Constraining the dark energy with galaxy clusters X-ray data 
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The equation of state characterizing the dark energy component is constrained by combining 
Chandra observations of the X-ray luminosity of galaxy clusters with independent measurements 
of the baryonic matter density and the latest measurements of the Hubble parameter as given by 
the HST key project. By assuming a spatially flat scenario driven by a "quintessence" component 
with an equation of state p x — uop x we place the following limits on the cosmological parameters 



u and il m : (i) —1 < uj < —0.55 and Q B 



0.32 



+0.027 



(la) if the equation of state of the dark 



energy is restricted to the interval —1 < uj < (usual quintessence) and (ii) uj — — 1.29to'792 
and Q m = 0.31±°;°^ (la) if uj violates the null energy condition and assume values < —1 (ex- 
tended quintessence or "phantom" energy). These results are in good agreement with independent 
studies based on supernovae observations, large-scale structure and the anisotropies of the cosmic 
background radiation. 

PACS numbers: 98.80-k; 98.80.Es; 98.65.Cw 



I. INTRODUCTION 

Recent observations of Type la supernovae (SNe la) 
have provided direct evidence that the universe may be 
accelerating Q. These results, when combined with 
measurements of cosmic microwave background radia- 
tion (CMB) anisotropies and dynamical estimates of 
the quantity of matter in the Universe, suggest a spa- 
tially flat universe composed of nearly ~ 1/3 of mat- 
ter (baryonic + dark) and ~ 2/3 of an exotic compo- 
nent endowed with large negative pressure, the so-called 
"quintessence". Nowadays, it has been promptly recog- 
nized that the question related to the nature of this dark 
energy is one of the most challenging problems of modern 
astrophysics, cosmology and particle physics. 

The absence of a natural guidance from particle physics 
theory about the nature of this dark component gave 
origin to an intense debate and many theoretical spec- 
ulations. In particular, a cosmological constant (A) - 
the most natural candidate - is the simplest but not the 
unique possibility; A is a time independent and spatially 
uniform dark component, which is described by a per- 
fect fluid with p v = —p v . Some other candidates ap- 
pearing in the literature are: a decaying vacuum energy 
density, or a time varying A-term 0], a time varying 
relic scalar field component (SFC) which is slowly rolling 
down its potencial the so-called "X-matter", an ex- 
tra component simply characterized by an equation of 
state p x = ujp x (XCDM) 0,0, the Chaplygin gas whose 
equation of state is given by p — —A/p where A is a 
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positive constant [6|, models based on the framework of 
brane- induced gravity , among others Q • For SFC and 
XCDM scenarios, the to parameter may be a function of 
the redshift (see, for example, |3) or still, as has been re- 
cently discussed, it may violate the null energy condition 
and assume values < — 1 |l(J . 

In order to improve our understanding of the actual 
nature of the dark energy, an important task nowadays 
in cosmology is to find new methods or to revive old ones 
that could directly or indirectly quantify the amount of 
dark energy present in the Universe, as well as determine 
its equation of state. In this concern, the possibility of 
constraining cosmological parameters from X-ray lumi- 
nosity of galaxy clusters constitutes an important and 
interest ing tool. This method was originally proposed by 
Sasaki [ll| and Pen [l2 with basis on measurements of 
the mean baryonic mass fraction in clusters as a func- 
tion of redshift. A recent application of a new version 
of this test was performed by Allen et al. 0, ^| (see 
also 0|) wn0 analyzed the X-ray observations in some 
relaxed lensing clusters observed with Chandra spanning 
the redshift range 0.1 < z < 0.5. By inferring the cor- 
responding gas mass fraction they placed observational 
limits on the total matter density parameter, f2 m , as well 
as on the density parameter, associated to the vac- 
uum energy density. 

In the present paper, by following the methodology 
presented in |14|. we discuss quantitatively how the ob- 
servations of X-ray gas mass fraction of galaxy clusters 
constrains the cosmic equation of state describing the 
dark energy component. In order to detect the pos- 
sibility of bias in the parameter determination due to 
the imposition uj > — 1 we have studied two different 
cases: the usual quintessence (—1 < to < 0) and the ex- 
tended quintessence (also named "phantom" energy 10']), 
in which the 10 parameter may assume values < — 1. 
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In the former case, a good agreement between theory 
and observations is possible if 0.3 < f2 m < 0.35 (68% 
c.l.) and ijj < —0.55. These results are in line with re- 
cent analyzes from distant SNe la [ljj, SNe + CMB @, 
gravitational lensing statistics ^} an d the existence of 
old high redshift objects (OHRO's) 0. For extended 
quintessence we obtain —2.1 < uj < —0.6 (68% c.l.) 
with the matter density parameter ranging in the interval 
0.27 < tt m < 0.34 (68% c.l.). 

This paper is organized as follows. In Section II we 
present the basic field equations and the distance for- 
mulas relevant for our analysis. The corresponding con- 
straints on the cosmological parameters to and f2 m are 
investigated in Section III. We finish the paper by sum- 
marizing the main results in the conclusion Section. 



II. BASIC EQUATIONS 

For a spatially flat, homogeneous, and isotropic cos- 
mologies driven by nonrelativistic matter and a sep- 
arately conserved exotic fluid with equation of state, 
p x = u)px, the Friedman's equation is given by: 



(|) 2 = Hi 



,R > 



,R \ 



n m (^) B + (i-fim)(^) 3(1+w) 
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FIG. 1: The model function /™° d (Eq. 8) as a function of the 
redshift for selected values of ui and fixed values of fi m = 0.3, 
Q, h h 2 = 0.0205 and h = 0.72. 



where an overdot denotes derivative with respect to time 
and H a = 100/iKm.s _1 .Mpc _1 is the present value of the 
Hubble parameter. 

In order to derive the constraints from X-ray gas mass 
fraction in the next Section we shall use the concept of 
angular diameter distance, Da(z). Such a quantity can 
be easily obtained in the following way: consider that 
photons are emitted by a source with coordinate r = r\ 
at time t\ and are received at time t a by an observer 
located at coordinate r = 0. The emitted radiation will 
follow null geodesies so that the comoving distance of the 
source is defined by (c = 1) 



which provides, when combined with Eq. (3), 



n DE 



H- 



dx 



(i + z) Ji x^/n^x- 1 + (i - o m )x-( 1 + 3 -) ' 



(5) 

For the standard cold dark matter model (SCDM) we 
set lu = in Eq. (3) and the angular diameter distance 
reduces to 



n scDM 



2H- 



(1 + z)3/2 



(l + z) 1/2 -l 



(6) 



ri 



dt 



ti 



dR 



R{t) J R (t) R{t)R(t) ' 



(2) 



By considering the above equations, it is straightforward 
to show that the comoving distance ri(z) can be written 
as 



n(z) 



1 



H„R„ 



dx 



^V^- 1 + (i - a„)x-( i + 3 -) ' 



(3) 

where the subscript o denotes present day quantities and 
x = Spp- = (l + z) -1 is a convenient integration variable. 
The angular diameter distance to a light source at r = t\ 
and t = t% which is observed at r = and t — t Q is 
defined as the ratio of the source diameter to its angular 
diameter, i.e., 



D A = - = R(ti)n = (1 + z)- 1 i? Q r(z), 



(4) 



III. CONSTRAINTS FROM X-RAY GAS MASS 
FRACTION 

In our analysis we consider the Chandra data analysed 
in recent papers by Allen et al. 0, and Schmidt et 
ai.jl^. The specific data set consists of six clusters dis- 
tributed over a wide range of redshifts (0.1 < z < 0.5). 
The clusters studied are all regular, relatively relaxed sys- 
tems for which independent confirmation of the matter 
density parameter results is available from gravitational 
lensing studies. As discussed in Ref. 0], the system- 
atic uncertainties are < 10% (i.e. typically smaller than 
the statistical uncertainties). The X-ray gas mass frac- 
tion (/gas) values were determined for a canonical radius 
r 2500i which is defined as the radius within which the 
mean mass density is 2500 times the critical density of 
the Universe at the redshift of the cluster. Two data sets 
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FIG. 2: Confidence regions in the fi m — to plane by assuming 
the SCDM model as the default cosmology. The regions in the 
graph correspond to 68%, 95% and 99% likelihood contours 
for flat quintessence scenarios. 



were generated from these data. One in which the SCDM 
model with H = 50Km.s _1 .Mpc _1 is used as the default 
cosmology and the other one in which the default cosmol- 
ogy is the ACDM scenario with H a — 70Km.s _1 .Mpc _1 , 
f2 m = 0.3 and f^A = 0.7. In what follows we constrain the 
basic cosmological parameters using the SCDM scenario 
as the default cosmology. 

By assuming that the baryonic mass fraction in galaxy 
clusters provides a fair sample of the distribution of 
baryons at large scale, the matter content of the universe 
can be expressed as [2(1 |21j : 



/ gas (l + 0.19^/2) : 



(7) 



where stands for the baryonic mass density parame- 



ter. Since / gas oc D A 
by 



3/2 



111! , the model function is defined 



/gas (^i) 



0.19/1 3 / 2 ) n n 



2h 



D s A CDM ( Zi ) 



1.5 



(8) 



where the term (2ft,) 3 / 2 represents the change in 
the Hubble parameter between the default cos- 
mology and quintessence scenarios while the ratio 
D A CBM (zi)/D A )E (zi) accounts for deviations in the ge- 
ometry of the universe from the SCDM model. Figure 1 
shows the behavior of /™° d as a function of the redshift 
for some selected values of co having the values of fib and 
h fixed. The value of O m is fixed at 0.3 as suggested 



FIG. 3: The same as in Fig. 2 by assuming the Gaussian 
priors h = 0.72 ± 0.08 and Q h h 2 = 0.0205 ± 0.0018. 



by dynamical estimates on scales up to about 2/i _1 Mpc 
|22j . For the sake of comparison, the current favored 
cosmological model, namely, a flat scenario with 70% of 
the critical energy density dominated by a cosmological 
constant (ACDM) is also shown. 

In order to determine the cosmological parameters il m 
and co we use a \ 2 minimization for the range of f2 m and 
co spanning the interval [0,1] in steps of 0.02 



i=l 



/gas, i] 



(9) 



where Uf asi are the symmetric root-mean-square errors 
for the SCDM data. The 68.3% and 95.4% confidence 
levels are defined by the conventional two-parameters x 2 
levels 2.30 and 6.17, respectively. 

In Fig. 2, by fixing the values of f2b (0.0205) and h 
(0.72), we show contours of constant likelihood (95% 
and 68%) in the £l m -cu plane. Note that the allowed 
range for both f2 m and co is reasonably large, showing 
the impossibility of placing restrictive limits on these 
quintessence scenarios from the considered X-ray gas 
mass fraction data. The best-fit model for these data 
occurs for f2 m = 0.33 and co = —1.0 with \ 2 = 1-98. 
Such limits become slightly more restrictive if we as- 
sume some a priori knowledge of the value of the product 
fl h h 2 = 0.0205 ± 0.0018 and of the value of the Hub- 
ble parameter h = 0.72 ± 0.08 0. To illustrate these 
new results, in Fig. 3 we show the confidence regions in 
the Q, m -co plane by assuming such priors. In this case, 
the best-fit model occurs for f2 m = 0.32, co = — 1 and 
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TABLE I: 


Limits to fi m 


and lo 




Method 


Reference 




LO 


V^lvlJj -f- DIN 6 la 


14] 


— U.o 


— — u.o 




rnl 

I9| 




< — 0.6 


blNe la 


|lb| 




< — U.00 


blNc la + \jL 


\Sb\ 


U.24 


< —0.7 


GL 


[17] 




—0.55 


OAT« T« i TOG 

blNe la + Lob 


I37| 




< —O.D 


Various 


I29| 


U.z — U.o 


< —O.D 


OTTRD'c 




u . o 


n 97 


CMB 


[381 


0.3 


< -0.5 




[32] 




< -0.96 


A6 


[331 


0.2-0.4 


< -0.5 


*(*) 


[301 


0.2 


~ -1.0 


CMB + SNe + LSS. 


[341 


0.3 


< -0.85 


CMB + SNe + LSS. 


[281 




< -0.71 


CMB + SNe + LSS." 


[281 




> -2.68 


SNe Ia a 


[281 


0.45 


-1.9 


SNe + X-ray Clusters a . 


[351 


~ 0.29 


-0.95 


X-ray Clusters 


This Paper 


~ 0.32 


< -0.5 


X-ray Clusters a 


This Paper 


~ 0.31 


-1.29 



"extended quintessence 



FIG. 4: Constraints on the fi m — lo plane for extended 
quintessence. The regions in the graph correspond to 68%, 
95% and 99% confidence limits. As in Fig. 3, Gaussian priors 
on the values of Q b h 2 and h were assumed. 



Xmin = 1-95 with the 1-a limits on 



and S7 m given, re- 



spectively, by 



and 



lu < -0.555, 



On = 0.320 



+0.027 
-0.014- 



For the sake of completeness, we also verified that by 
fixing lo = — 1 and extending the analysis for arbitrary 
geometries the results of 0] are fully recovered. 

So far we have assumed that the dark energy equation 
of state is constrained to be w > —1. However, as has 
been observed recently a dark component with lo < — 1 
appears to provide a better fit to SNe la observations 
than do ACDM scenarios (lo = —1) [HJ- I n f ac t; al- 
though having some unusual properties, this "phantom'' 
behavior is predicted by several scenarios as, for exam- 
ple, kinetically driven models [^J an d some versions of 
brane world cosmologies [2|| (see also [22} and references 
therein). In this concern, a natural question at this point 
is: how does this extension of the parameter space to 
lo < — 1 modify the previous results? To answer this 
question in Fig. 4 we show the 68% and 95% confidence 
regions in the "extended" Q m — lo plane by assuming the 
same a priori knowledge of the product fl^h 2 and of the 
value of the Hubble parameter as done earlier. From 
this analysis, we find O m = 0.312^°°^, lo = -1.29^;^ 



and x r 2 n in = 1-77 both results at 1-a level. By assum- 
ing no a priori knowledge on fl^h 2 and h we obtain 
to = —1.28^0 809 wn il e th e value of fi m remains approxi- 
mately the same. These limits should be compared with 
the ones obtained by Hannestad & Mortsell |23 by com- 
bining CMB + Large Scale Structure (LSS) + SNe la 
data. At 95.4% c.l. they found -2.68 < lo < -0.78. 

At this point we compare our results with other recent 
determinations of lo derived from independent methods. 
For example, for the usual quintessence (i.e., lo > — 1), 
Garnavich et al. 0] used the SNe la data from the 
High-z Supernova Search Team to find lo < —0.55 (95% 
c.l.) for flat models whatever the value of f2 m whereas for 
arbitrary geometry they obtained lo < —0.6 (95% c.l.). 
Such results agree with the constraints obtained from a 
wide variety of different phenomena, using the "concor- 
dance cosmic" method |2||. In this case, the combined 
maximum likelihood analysis suggests lo < —0.6, which 
rules out an unknown component like topological defects 
(domain walls and string) for which lo = — ^ , being n the 
dimension of the defect. Recently, Lima and Alcaniz pffll ] 
investigated the angular size - redshift diagram (6(z)) 
in quintessence models by using the Gurvits' published 
data set [Ul- Their analysis sugests —1 < lo < —0.5 
whereas Corasaniti and Copeland found, by using 
SNe la data and measurements of the position of the 
acoustic peaks in the CMB spectrum, — 1 < lo < —0.93 
at 2<7. More recently, Jain et al. [33 used image sep- 
aration distribution function (A9) of lensed quasars to 
obtain —0.8 < lo < —0.4, for the observed range of 
Q m ~ 0.2 — 0.4 while Chae et al. 17] used gravitational 
lens (GL) statistics based on the final Cosmic Lens All 



Sky Survey (CLASS) data to find lo < -0.55] 



(68% 
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c.L). Bean and Melchiorri [34| obtained u> < —0.85 from 
CMB + SNe la + LSS data, which provides no significant 
evidence for quintessential behaviour different from that 
of a cosmological constant. A similar conclusion was also 
obtained by Schuecker et al. |35| from an analysis involv- 
ing the REFLEX X-ray cluster and SNe la data in which 
the condition uj > — 1 was relaxed. A more extensive list 
of recent determinations of the quintessence parameter oj 
is presented in Table I. 

IV. CONCLUSION 

The determination of cosmological parameters is a cen- 
tral goal of modern cosmology. We live in a special mo- 
ment where the emergence of a new "standard cosmol- 
ogy" driven by some form of dark energy seems to be 
inevitable. The uncomfortable situation for some comes 
from the fact that the emerging model is somewhat more 
complicated physically speaking, while for others it is 
exciting because although preserving some aspects of the 
basic scenario a new invisible actor which has not been 
predicted by particle physics is coming into play. 

Using the reasonable ansatz of constant gas mass frac- 
tion at large scale, we placed new limits on the f2 m and 
uj parameters for a flat dark energy model. The galaxy 
cluster data used corresponds to regular, relaxed systems 
whose /gas(f) profiles are essentially flat around r25oo> the 
mass results were confirmed from gravitational lensing 
studies and the residual systematic uncertainties in the 
/gas values are small ^4|. Naturally, the analysis pre- 
sented here also reinforces the interest in searching for 



X-ray data both for less relaxed clusters, and perhaps 
more important, at higher redshifts. Hopefully, our con- 
straints will be more stringent when further Chandra, 
XMM-Newton and more accurate gravitational lensing 
data for clusters become available near future. In this 
concern, we recall that X-ray data from galaxy clusters 
at high redshifts and the corresponding constraints for 
£l m will play a key hole in the coming years because 
their relative abundance (and consequently the value of 
O m itself) may also independently be checked trough the 
Sunyaev-Zeldovich effect [39j|. 

As we have seen, the X-ray data at present also favor 
eternal expansion as the fate of the Universe in accor- 
dance with SNe la data.!]. Our estimates of f2 m and 
uj are compatible with the results obtained from many 
independent methods (see Table I). We emphasize that 
a combination of these X-ray data with different meth- 
ods is very welcome not only because of the gain in pre- 
cision but also because most of cosmological tests are 
endowed with a high degree of degeneracy and may con- 
strain rather well only specific combinations of cosmolog- 
ical parameters but not each parameter individually. The 
basic results combining different methods will appear in 
a forthcoming communication |40|. 
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